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This short review is an introduction to a great variety of methods, the 
collection of which is called the Painleve analysis, intended at producing all 
kinds of exact (as opposed to perturbative) results on nonlinear equations, 
whether ordinary, partial, or discrete. 

1 The classical program of the Painleve school 
and its achievements 

It is impossible to understand anything to the Painleve property without keep- 
ing in mind the original problem as stated by L. Fuchs, Poincare and Painleve: 
to define new functions from ordinary differential equations (ODEs). This 
simply formulated problem implies to select those ODEs whose general solu- 
tion can be made singlevalued by some uniformization procedure (cuts, Rie- 
mann surface), so as to fit the definition of a function. This property (the 
possibility to uniformize the general solution of an ODE), nowadays called 
the Painleve property (PP), is equivalent to the more practical definition. 

Definition 1 The Painleve property of an ODE is the absence of movable 
critical singularities in its general solution. 

Let us recall that a singularity is said movable (as opposed to fixed) if its 
location depends on the initial conditions, and critical if multivaluedness takes 
place around it. Other definitions of the PP, excluding for instance the es- 
sential singularities, or replacing "movable critical singularities" by "movable 
singularities other than poles" , or "its general solution" by "all its solutions" , 
are incorrect. Two examples taken from Chazy jll] explain why this is so. 
The first example is the celebrated Chazy's equation of class III 

u'" -2uu" + 3u' 2 = 0, (I) 

whose general solution is only defined inside or outside a circle characterized 
by the three initial conditions (two for the center, one for the radius); this 
solution is holomorphic in its domain of definition and cannot be analytically 
continued beyond it. This equation therefore has the PP, and the only singu- 
larity is a movable analytic essential singular line which is a natural boundary. 
The second example |TT|, p. 360] is the third order second degree ODE 

(u'" - 2u'u") 2 + 4u" 2 (u" - u' 2 - 1 ) = 0, (2) 

whose general solution is singlevalued, 

u = e c l x+c 2/ci + ^A x + c ^ (3) 

but which also admits a singular solution (envelope solution) with a movable 
critical singularity, 

u = C2 — Logcos(a; — C\). (4) 
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For more details, see the arguments of Painleve |0, §2.6] and Chazy [jl?], §5.1]. 

The PP is invariant under an arbitrary homography on the dependent 
variable and an arbitrary change of the independent variable {nomographic 
group) 



(it, x) i— > (U,X), u(x) 



a(x)U(X)+/3(x) 



j(x)U(X) + S{x) ' 
(a, 0, 7, S, £) functions, ad — f3j ^ 0. 



x = ax), 



(5) 



Every linear ODE possesses the PP since its general solution depends lin- 
early on the movable constants, so, in order to define new functions, one must 
turn to nonlinear ODEs in a systematic way: first order algebraic equations, 
then second order, . . . The current achievements are the following. 

First order algebraic ODEs (polynomial in u, v! , analytic in x) define only 
one function, the Weierstrass elliptic function p, new in the sense that its 
ODE 

u' 2 — 4u 3 + §2U + 53 = 0, (32,53) arbitrary complex constants, (6) 

is not reducible to a linear ODE. Its only singularities are movable double 
poles. 

Second order algebraic ODEs (polynomial in u,u',u", analytic in x) de- 
fine six functions, the Painleve functions Pn, n — 1, • ■ • , 6, new because they 
are not reducible to either a linear ODE or a first order ODE. This question 
of irreducibility, the subject of a long dispute between Painleve and Joseph 
Liouville, has been rigorously settled only recently |)l| . The canonical repre- 
sentatives of P1-P6 in their equivalence class under the group (pi) are 



PI 
P2 

P3 
P4 
P5 
P6 



u 
u" 



Qu z - 
2m 3 - 

u 

2u ' 
1 

2m 



x, 

xu - 

r/ 

x 
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2 U 



au 



7U 



Ax 2 



fi_ 

4;r; 



Axu z 



2x 2 u ■ 



S 

2au + — i 
u 



1 



u = — 



u-1 
1 



vf_ (u- l) 2 

X x z 



au 



u — 1 u 
u(u — l)(u — x) 
x 2 (x-l) 2 



/2 

u — 



1 



1 

x x — 1 
x - 1 



u r u(u + 1) 
x u — 1 



1 



'7 



u — X 

5 x(x - 1) 

(u — l) 2 (u — x) 2 



in which a, j3, 7, S are arbitrary complex parameters. Their only singularities 
are movable poles (in the e x complex plane for P3 and P5, in the x plane for 
the others), with in addition three fixed critical singularities for P6, located 
at x = 00, 0, 1. 

Third and higher order ODEs jy], 0, |35|, |3(| have not yet defined new 
functions. Although there are some good candidates (the Gamier system 
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|47|1 , several fourth order ODEs |36| g5|, which all have a transcendental de- 
pendence on the constants of integration) , the question of their irreducibility 
(to a linear, Weierstrass, or Painleve equation) is very difficult and still open. 
To understand the difficulty, it is sufficient to consider the fourth order ODE 
for u(x) defined by 

u = u\ + U2, u'{ — 6uf + x, u 2 ' = 61*2 + x. (7) 

This ODE (easy to write by elimination of u\,U2) has a general solution 
which depends transcendentally on the four constants of integration, and it is 
reducible. 

The master equation P6 was first written by Picard in 1889 in a particular 
case, in a very elegant way. Let ip be the elliptic function defined by 

Lp : y^Lp(y,x), y = f - ^ ==, (8) 

Joo \JZ{Z - 1)(Z - X) 

and let loi(x),uj2{x) be its two half-periods. Then the function 

u : x 1-^ u(x) — Lp(2ciUJi(x) + 2c2UJ2{x),x), (9) 

with (01,02) arbitrary constants, has no movable critical singularities, and it 
satisfies a second order ODE which is P6 in the particular case a = (3 = 7 = 
5 — 1/2 = 0. The generic P6 was found simultaneously from two different ap- 
proaches, the nonlinear one of the Painleve school as said above |8~i]| , and the 
linear one of R. Fuchs [Q as an isomonodromy condition. In the latter, one 
considers a second order linear ODE for tp(t) with four Fuchsian singularities 
of crossratio x (located for instance at t = 00, 0, 1, x), with in addition, as pre- 
scribed by Poincare for the isomonodromy problem, one apparent singularity 
t = u, 

2d 2 ijj A B C D 3 



ip At 2 i 2 {t-l) 2 {t-xf t{t-l) A{t- 

a b 



' t(t-l){t-x) ' t(t-l)(t-u)' ^ 

(A, B,C, D) denoting constants and (a, b) parameters. The requirement that 
the monodromy matrix (which transforms two independent solutions "01 , V'2 
when t goes around a singularity) be independent of the nonapparent singu- 
larity x results in the condition that it, as a function of x, satisfies P6. 

A useful by-product of this search for new functions is the construction 
of several exhaustive lists {classifications) of second |^. [|, |3l], [3^|, third 
(ll], ^, p6| |, fourth [j?], pBq ), or higher order ODEs, whose general solution 
is explicitly given because they have the PP. Accordingly, if one has an ODE 
in such an already studied class (e.g. second order second degree binomial- 
type ODEs |3lJ u" 2 ~ F(u',u,x) with F rational in u' and u, analytic in x), 
and which is suspected to have the PP (for instance because one has been 
unable to detect any movable critical singularity, see Section |^), then two 
cases are possible: either there exists a transformation (^|) mapping it to a 
listed equation, in which case the ODE has the PP and is explicitly integrated, 
or such a transformation does not exist, and the ODE has not the PP. 
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2 Integrability and Painleve property for par- 
tial differential equations 

Defining the PP for PDEs is not easy, but this must be done for future use in 
sections || (the Painleve test) and || (proving the PP). Such a definition must 
involve a global, constructive property, which excludes the concept of general 
solution. Indeed, this is only in nongeneric cases like the Liouville equation 
that the general solution of a PDE can be built explicitly. This is where the 
Backhand transformation comes in. Let us first recall the definition of this 
powerful tool (for simplicity, but this is not a restriction, we give the basic 
definitions for a PDE defined as a single scalar equation for one dependent 
variable u and two independent variables [x,i)). 

Definition 2 (Refs. vol. Ill chap. XII], J&jj/j A Backlund transfor- 
mation (BT) between two given PDEs 

E x {u,x,t) =0, E 2 (U,X,T) = Q (11) 

is a pair of relations 

Fj(u,x,t,U,X,T) = 0, j = 1,2 (12) 

with some transformation between (x,t) and (X,T), in which Fj depends on 
the derivatives ofu(x,t) andU(X,T), such that the elimination of u (resp. U) 
between (F\,F<i) implies E2(U, X,T) = (resp. Ei(u,x,t) — 0). In case the 
two PDEs are the same, the BT is also called the auto-BT. 

Under a reduction PDE^ODE, the BT reduces to a birational transfor- 
mation (also initials BT!), which is not involved in the definition of the PP for 
ODEs. Therefore one needs an intermediate (and quite important) definition 
before defining the PP. 

Definition 3 A PDE in N independent variables is integrable if at least 
one of the following properties holds. 

1. Its general solution can be obtained, and it is an explicit closed form 
expression, possibly presenting movable critical singularities. 

2. It is linearizable. 

3. For N > 1, it possesses an auto-BT which, if N — 2, depends on an 
arbitrary complex constant, the Backlund parameter. 

4-. It possesses a BT to another integrable PDE. 

Examples of these various situations are, respectively: the PDE u x u t + 
uu xt = with general solution u = y/ f(x) + g(t), which presents movable 
critical singularities and can be transformed into the d'Alembert equation; 
the Burgers PDE u t + u xx + 2uu x — 0, linearizable into the heat equation 
ipt + ^xx = 0; the KdV PDE u t + u xxx — 6uu x = 0, which is integrable by 
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the inverse spectral transform (1ST); the Liouville PDE u x t + e u = 0, which 
possesses a BT to the d'Alembert equation tp xt = 0. 

We now have enough elements to give a definition of the Painleve property 
for PDEs which is indeed an extrapolation of the one for ODEs. 

Definition 4 The Painleve property (PP) of a PDE is its integrability (defi- 
nition [|) and the absence of movable critical singularities near any nonchar- 
acteristic manifold. 

One will retain that the Painleve property is a more demanding property 
than the mere integrability. 

The PP for PDEs is invariant under the natural extension of the homo- 
graphic group (|5|), and classifications similar to those of ODEs have also been 
performed for PDEs, in particular second order first degree PDEs |33[ [j4| . 
isolating only the already known PDEs (Burgers, Liouville, sine-Gordon, 
Tzitzeica, etc). Classifications based on other criteria, such as the existence 
of an infinite number of conservation laws pq], isolate more PDEs, which are 
likely all integrable in the sense of definition^; it would be interesting to check 
that, under the group of transformations generated by Backlund transforma- 
tions and hodograph transformations, each of them is equivalent to a PDE 
with the PP. 

If one performs a hodograph transformation (typically an exchange of the 
dependent and independent variables u and x) on a PDE having the PP, the 
transformed PDE possesses a weaker form of the PP in which, for instance, all 
leading powers and Fuchs indices become rational numbers instead of integers. 
Details can be found e.g. in Ref. H| . For instance, the Harry-Dym, Camassa- 
Holm and DHH equations |39|| can all be mapped to a PDE with the PP 
by some hodograph transformation. 

If the above definition of the PP for PDEs is really an extrapolation of 
the one for ODEs, then, given a PDE with the PP, every reduction to an 
ODE which preserves the differential order (i.e. a noncharacteristic reduction) 
yields an ODE which necessarily has the PP. This proves the conjecture of 
Ablowitz, Ramani and Segur provided the definition |] really extrapolates 
that for ODEs, and this is the difficult part of this question. 

There are plenty of such examples of reductions, for instance the self-dual 
Yang-Mills equations admit reductions to all six Pn equations [ |64[ . 

Given a PDE, the process to prove its integrability is twofold. One must 
first check whether it may be integrable, for instance by applying the Painleve 
test, as will now be explained in section |^. Then, in case of a nonnegative 
answer, one must build explicitly the elements which are required to establish 
integrability, for instance by using methods described in section 

Although partially integrable and nonintegrable equations, i.e. the major- 
ity of physical equations, admit no BT, they retain part of the properties of 
(fully) integrable PDEs, and this is why the methods presented below apply 
to both cas es as well. One such example is included below for information, 
section ff.5[ 
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3 The Painleve test for ODEs and PDEs 

The test is an algorithm providing a set of necessary conditions for the equa- 
tion to possess the PP. Its full, detailed version can be found in Ref. Jl?], 
§6], and below is only given a short presentation of its main subset, known 
as the method of pole-like expansions, due to Kowalevski and Gambier fl45| . 
The generated necessary conditions are a priori not sufficient, this was proven 
by Picard (1894), who exhibited the example of the ODE with the general 
solution u = p(ALog(x — Ci) + 02,32,53), namely 

u" — —3 — — i^- 9 })- - U ' 2 = 0, (13) 

Au 6 - g 2 u - g 3 V 2 J X^/Au 3 - g 2 u - g 3 

which has the PP iff 2niX is a period of the elliptic function p, a transcendental 
condition on (A, g 2l gz) impossible to obtain in a finite number of algebraic 
steps such as the Painleve test. Therefore it is wrong to issue statements like 
"The equation passes the test, therefore it has the PP" . The only way to prove 
the PP is, 

- for an ODE, either to explicitly integrate with the known functions (so- 

lutions of linear, elliptic, hyperelliptic (a generalization of elliptic) or 
Painleve equations), or, if one believes to have found a new function, 
to prove both the absence of movable critical singularities and the irre- 
ducibility [ j9l) to the known functions, 

- for a PDE, to explicitly build the integrability elements of the definition |J. 

Let us return to the test itself. It is sufficient to present the method of 
pole-like expansions for one iVth order equation 

E(x,t,u) = 0, (14) 

in one dependent variable u and two independent variables x, t. Movable 
singularities lay on a codimension one manifold 

<p(x,t)-ipo = 0, (15) 

in which the singular manifold variable tp is an arbitrary function of the in- 
dependent variables and tpo an arbitrary movable constant. Basically p4| , 
the test of Weiss, Tabor and Carnevale consists in checking the existence of 
all possible local representations, near ip(x, t) — ipo = 0, of the general solu- 
tion (whatever be its definition, difficult for PDEs) as a locally single valued 
expression, e.g. the Laurent series 

+oc +00 

« = J2 u ^ j+p ' -p e ^> E = E E ^ j+ ^ -1 e ^> ( 16 ) 

3=0 1=0 

with coefficients Uj, Ej independent of the expansion variable \- The natural 
choice x — V ~ ^0 [p4j generates lengthy expressions Uj , Ej . Fortunately, 
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there is much freedom in choosing only required to vanish as tp — po and to 
have a homographic dependence on tp — ipo so as not to alter the structure of 
movable singularities, hence the result of the test. The unique choice which 
minimizes the expressions and puts no constraint on ip is |16| 



<P ~ Po 

X = 



Pxx i \ 

fx - -X {f - fo) 

2fx 



-\ -1 

fx fxx 



P-PO ^fx 



fx + 0, (17) 



in which x denotes an independent variable whose component of grad p does 
not vanish. The expansion coefficients Uj,Ej are then invariant under the 
six-parameter group of homographic transformations 

<^!t^7> a'd'-b'c'^O, (18) 
dp + a' 

in which a',b',c',d' are arbitrary complex constants, and these coefficients 
only depend on the following elementary differential invariants and their 
derivatives: the Schwarzian 

a r i Pxxx 3 / p> xx \ 

S = {p;x} = — — , (19) 

Px 2 \Px J 

and one other invariant per independent variable t, y, . . . 

C = -tpt/Px, K = ~p y /p x , ... (20) 

The two invariants S, C are linked by the cross-derivative condition 

X = {{pxxx)t - (<Pt)xxx)/Px — S t + C xxx + 2C X S + CS X = 0, (21) 

identically satisfied in terms of p. 

For the practical computation of (uj,Ej) as functions of (S,C) only, 
i.e. what is called the invariant Painleve analysis, the variable tp disappears, 
and the only information required is the gradient of the expansion variable \i 

Xx = l + fx 2 , Xt = -C + C xX - \{CS + C xx )x 2 - (22) 

with the constraint ( |2l| ) between S and C. 

Consider for instance the Kolmogorov-Petrovskii-Piskunov (KPP) equa- 
tion g Q 

E(u) = bu t - u xx + juu x + 2gT 2 (« - ei)(u - e2){u - e 3 ) = 0, (23) 

with (b,"f,d 2 ) real and ej real and distinct, encountered in reaction-diffusion 
systems (the convection term uu x j87| is quite important in physical applica- 
tions to prey-predator models). 

The first step, to search for the families of movable singularities u ~ 
uqx p , E ~ E()X q , uq 7^ 0, results in the selection of the dominant terms E(u) 

E(u) = -u xx + -fuu x + 2d~ 2 u 3 , (24) 
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which provide two solutions (p, uq) 

p=-l,q=-3,-2-ju + 2d~ 2 ul = 0. (25) 

The necessary condition that all values of p be integer is satisfied. 

The second step is, for every selected family, to compute the linearized 
equation, 

_ E(u + ew) - E(u) 
(E (u))w — hm 

= {-d 2 x +^ud x + 1 u x + 6d- 2 u 2 )w = 0, (26) 

then its Fuchs indices i near x = as the roots of the indicial equation 

P(i) = lim x^i-dl + lu oX p d x + ipuox 13 ' 1 + M~ 2 ul X 2p )x t+p (27) 

= -{i-l){i-2)+ju {i-2) + 6d- 2 ul (28) 
= -(j + i)(i-4-7Uo) = 0, (29) 

and finally to enforce the necessary condition that, for each family, these two 
indices be distinct integers jl3[ |l^] . Considering each family separately would 
produce a countable number of solutions, which is incorrect. Considering the 
two families simultaneously, the diophantine condition that the two values 
i\, ii of the Fuchs index 4 + 71*0 be integer has a finite number of solutions, 
namely || App. I] 

j 2 d 2 = 0, (*!, ^ = (4, 4), u = (-d, d), (30) 
7 V = 2, (h,i 2 ) = (3, 6), 7 uo = (-1, 2), (31) 
7 2 d 2 = -18, (ii, i 2 ) = (-2, 1), 7U0 = (-6, -3). (32) 

It would be wrong at this stage to discard negative integer indices. Indeed, in 
linear ODEs such as ( p6| ) , the single valuedness required by the Painleve test 
restricts the Fuchs indices to integers, whatever be their sign. Let us proceed 
with the first case only, 7 = (the usual KPP equation). 
The recurrence relation for the next coefficients Uj, 

Vj > 1 : E j = P(j) Uj + Q 3 -({u, \ l<j}) = (33) 

depends linearly on Uj and nonlinearly on the previous coefficients it;. 

The third and last step is then to require, for any admissible family and 
any Fuchs index i, that the no-logarithm condition 

Vz G Z, P(i) = : Qi = (34) 

holds true. At index i = 4, the two conditions, one for each sign of d fll5| , 

Qi = C[(bdC + sx - 3ei)(bdC + st - 3e 2 ){bdC + si - 3e 3 ) 

~ 3b 2 d 3 (C t + CC X )] =0, si - ex + e 2 + e 3 , (35) 

are not identically satisfied, so the PDE fails the test. This ends the test. 
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If instead of the PDE ( |23| ) one considers its reduction u(x, t) = U(£),£ = 
x — ct to an ODE, then C = constant = c, and the two conditions Q4 = 
select the seven values c = and c 2 = (si — 3ek) 2 (bd)~ 2 , k = 1,2,3. For 
all these values, the necessary conditions are then sufficient since the general 
solution U(£) is singlevalued (equation numbered 8 in the list of Gambier 
reproduced in |j54|). 

It frequently happens that the Laurent series ( [Trf ) only represents a par- 
ticular solution, for instance because some Fuchs indices are negative integers, 
e.g. the fourth order ODE [g p. 79] 

u"" + 3uu" - Au' 2 = (36) 

which admits the family 

p = -2, u = -60, Fuchs indices (-3, -2, -1, 20). (37) 

The series (|l^) depends on two, not four, arbitrary constants, so two are 
missing and may contain multivaluedness. In such cases, one must perform a 
perturbation in order to represent the general solution and to test the missing 
part of the solution for multivaluedness. 

This perturbation [jl9| is close to the identity (for brevity, we skip the t 
variable) 

+00 +00 
x unchanged, u = ^ e n u {n) : E = ^ e n E^ = 0, (38) 

n=0 ra=0 

where, like for the a— method of Painleve, the small parameter e is not in the 
original equation. 

Then, the single equation ( |l4|) is equivalent to the infinite sequence 

n = E { ^=E{x,u { ^) =0, (39) 
Vn > 1 E^=E'(x,u^)u^ + i? (n) (x,M (0) , . . . ,u (n_1) ) = 0, (40) 

with i?W identically zero. From a basic theorem of Poincare (Tt], Theorem II, 
§5.3], necessary conditions for the PP are 



the general solution of (39) has no movable critical points, 



the general solution of (40) has no movable critical points, 



- for every n > 2 there exists a particular solution of (40) without movable 
critical points. 

Order zero is just the original equation ( |l4| ) for the unknown u^°\ so one 
takes for the already computed (particular) Laurent series (|l6|). 
Order n = 1 is identical to the linearized equation 

E (1) = E'(x,u {0) )u w =0, (41) 

and one must check the existence of N independent solutions locally 
singlevalued near X = 0, where TV is the order of (|l4[). 
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The two main implementations of this perturbation are the Fuchsian per- 
turbative metho d [fl9f and the nonFuchsian perturbative method |fr"jf . In the 
above example (p6|), both methods indeed detect multivaluedness, at pertur- 
bation order n = 7 for the first one, and n — 1 for the second one (details 
below). 



3.1 The Fuchsian perturbative method 

Adapted to the presence of negative integer indices in addition to the ever 
present value — 1 , this method [lj| generates additional no-log conditions 
(34). Denoting p the lowest integer Fuchs index, p < — 1, the Laurent series 
for 

+ oo 

u«=£u$V* (42) 

3=P 

represents a particular solution containing a number of arbitrary coefficients 
equal to the number of Fuchs indices, counting their multiplicity. If this 
number equals N, it represents the general solution of (|4l]). Two examples 
will illustrate the method jl7|, §5.7.3]. 
The equation 

u" + 4uu' + 2u 3 = (43) 

possesses the single family 

p = -1, E^ = u£Vo 0) - !) 2 = 0, indices (-1.0), (44) 

with the puzzling fact that should be at the same time equal to 1 according 
to the equation Eq ^ = 0, and arbitrary according to the index 0. The necessity 
to perform a perturbation arises from the multiple root of the equation for 
Uq°\ responsible for the insufficient number of arbitrary parameters in the 
zeroth order series u^ '. The application of the method provides 

u (o) = x -i 

(the series terminates), (45) 

E'(x.u^) = d 2 x + 4 X - 1 d x + 2 X - 2 , (46) 

— Uq u q arbitrary, (47) 

= x- 2 (xV 2 ))" + 2 U « 2 X - 3 = 0, (48) 
t*W = -24 1)2 X - 1 (Log X -l). (49) 

The movable logarithmic branch point is therefore detected in a systematic 
way at order n = 2 and index i = 0. This result was of course found long ago 
by the a-method III §13, p 221]. 

The equation (Bq) possesses the two families 

p = -2, 4 0) = -60, ind. (-3, -2, -1, 20), E = u"" + 3uu" - 4u' 2 , (50) 
p = —3, itg ^ arbitrary, indices (— 1, 0), E = 3uu" — 4m' 2 . (51) 
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The second family has a Laurent series (p : +00) which happens to termi- 



nate [19 



,(0) 



c(x — xq) 3 — 60(x — xq) 2 , (c, xq) arbitrary. (52) 



For this family, the Fuchsian perturbative method is then useless, because the 
two arbitrary coefficients corresponding to the two Fuchs indices are already 
present at zeroth order. 

The first family provides, at zeroth order, only a two-parameter expansion 
and, when one checks the existence of the perturbed solution 



+00 

y 

71=0 



+00 

3=0 



(53) 



one finds that coefficients u^ , U_\, u_%, u_\ can be chosen arbitrarily, and, 
at order n = 7, one finds two violations fill 

^.(7) _ (0) (l) 7 n ^,(7) _ (0) 2 (l) 6 (1) n fKA\ 
Q-l = u 20 u -3 = °> Q20 = U 20 u -3 u -2 = °> ( 54 ) 

implying the existence of a movable logarithmic branch point. 

3.2 The nonFuchsian perturbative method 

Whenever the number of indices is less than the differential order of the equa- 
tion, the Fuchsian perturbative method fails to build a representation of the 
general solution, thus possibly missing some no-log conditions. The missing 



solutions of the linearized equation (41) are then solutions of the nonFuchsian 
type near X = 0. 



In section 3.1, the fourth order equation ( |36| ) has been shown to fail the 
test after a computation practically untractable without a computer. Let us 
prove the same result without computation at all |7l| . The linearized equation 

E« =E'{ Xl u^)u^ = [d 4 x +3u^d 2 x -8u^d x + 3u^]u (1 ^0, (55) 

is known globally for the second family because the two-parameter solution 



(52) is closed form, therefore one can test all the singular points \ °f (p£ 
These are x — (nonFuchsian) and \ — 00 (Fuchsian), and the key to the 
method is the information obtainable from \ = 00 ■ Let us first lower by 
two units the order of the linearized equation (|55|), by "subtracting" the two 
global single valued solutions — d Xa u^ and d c u^\ i.e. — x 4 ,X~ 3 i 

= X - 4 V : [dl - lSx^dx +3cx -3 ~ 60x~ 2 ]v" = 0, (56) 

Then the local study of % = 00 is unnecessary, since one recognizes the Bessel 
equation. The two other solutions in global form are 



c^0: v'i = X - 3 oF 1 (24;-3c/x) = X 17/2 J 23 ( V^c/x), (57) 
v'{ = X 17/2 N 23 (^12^), (58) 

where the hypergeometric function o-F\(24; — 3c/ x) is single valued and pos- 
sesses an isolated essential singularity at x = 0, while the fonction ^23 of 
Neumann is multivalued because of a Log x term. 
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4 Singularity-based methods towards integra- 
bility 

In this section, we review a variety of singularity-based methods able to pro- 
vide some global elements of integrability. The singular manifold method of 
Weiss et al. [Q is the most important of them, but it is not the only one. 
A prerequisite notion is the singular part operator T>, 

Log <p i— > V Log <p = ut{0) — u T (oo), (59) 

in which the notation «t(^o), which emphasizes the dependence on ipo, stands 
for the principal part of the Laurent series (|l6[), 

u T M = J £u jX j+v . (60) 

3=0 

In our KPP example fl25| ) with 7 = 0, this operator is T> = dd x . 



4.1 Linearizable equations 

When a nonlinear equation can be linearized, the singular part operator de- 
fined in (|5^) directly defines the linearizing transformation. 

For instance, the Kundu-Eckhaus PDE for the complex field U(x, t) |^| 



2 



iU t + aU xx + {^\U\ l + 2be^(\U\%)U = 0, (a, p, b, 7) € U, (61) 
a 

with a (3b cos 7 7^ 0, passes the test iff 21 b 2 = (3 2 . Under the parametric 
representation 

U = ^T x e l \ (62) 
the equivalent fourth order PDE for u |2l| 

a, 2 3 fl 2 - (6sin7) 2 4 

~^\U j xxxx'U j x 4~ ^xx 2v JX U xx U xxx ) 1 2 Vj x U xx 

2 a 

+2(b cos j)ulu xxx + — (uttul + u xx u\ - 2u t u x u xt ) = (63) 
2a 

admits two families, namely in the case b 2 = [3 2 , 

u ~ — — Log-0, indices —1,0,1,2, (64) 

2p cos 7 

3 

u^— Log-0, indices —3,-1,0,2, (65) 

2p cos 7 

in which (Log ip) x is the x of the invariant Painleve analysis. When the test 
is satisfied (b 2 = f3 2 ), the linearizing transformation [^9| is provided by |2lf| 
the singular part operator of the first family, which maps the nonlinear PDE 
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( JQlT) to the linear Schrodinger equation for V obtained by setting b = (3 = 
Kundu-Eckhaus([/) ^=^> Schr6dinger(F), iVt + aV xx — 0. (66) 



2/5 cos 7 

4.2 Auto-Backlund transformation of a PDE: the singu- 
lar manifold method 

Widely known as the singular manifold method or truncation method because 
it selects the beginning of a Laurent series and discards ("truncates") the 
remaining infinite part, this method was introduced by Weiss et al. J94| and 
later improved in many directions [ll , 46 ]7(], ^2], |2| [7^] . Its most recent 
version can be found in the lecture notes of a CIME school |l^, |68| , to which 
we refer for further details. 

The goal is to find the Backhand transformation or, if a BT does not 
exist, to generate some exact solutions. Since the BT is itself the result of an 
elimination [jl2| between the Lax pair and the Darboux involution, the task 
splits into the two simpler tasks of deriving these two elements. Let us take 
one example. 

The modified Korteweg-de Vries equation (mKdV) 

mKdV(w) = bw t + (w xx - 2w 3 /a 2 ) x = 0, (68) 

is equivalently written in its potential form 

p-mKdV(r) = br t + r xxx - 2r 3 Ja 2 + F(t) =0, w = r x , (69) 

a feature which will shorten the expressions to come. This last PDE admits 
the two opposite families (a is any square root of a 2 ) 

p = 0~, q = — 3, r~aLogV>, indices (—1, 0, 4), T> — a, (70) 

and the results to be found are: 

- the Darboux involution 

r = T>LogY + R, (71) 

a relation expressing the difference of two solutions r and R of p-mKdV 
as the logarithmic derivative VLogY, in which T> = a is the singu- 
lar part operator of either family, and Y is a Riccati pseudopotential 
equivalent to the Lax pair (see next item), 

- the Lax pair, here written in its equivalent Riccati representation, 

^=A(i-y)-2^, (72) 

y y a 

ttt.i (-^ + (2^ + 2^-4^) , (73) 
y y \ a a z a J x 

in which W satisfies the mKdV equation (^8[) and A is the spectral 
parameter, 
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the BT, by some elimination between the above two items. 



This program is achieved by defining the truncation 82 



r = aLogY + R, (74) 

in which r satisfies p-mKdV, R is a yet unconstrained field, Y is the most 
general homographic transform of x which vanishes as x vanishes, 

Y- 1 =B( X - 1 +A), (75) 



A and B are two adjustable fields, and the gradient of x is (22). The l.h.s. of 
the PDE is then 

6 

p-mKdV(r) = ^ Ej (S, C, A, B, R)Y j - 3 , (76) 

3=0 

and the system of determining equations to be solved is 

Vj E j (S,C,A,B,R) = 0. (77) 

This choice of Y (|7|) is necessary to implement the two opposite families 
feature of mKdV. The general solution of the determining equations intro- 
duces an arbitrary complex constant A and a new field W J82{ 

W = (R-aLogB) x , A= W/a, 
bC = 2W x /a- 2W 2 /a 2 +4A 2 , 
S = 2W x /a- 2W 2 /a 2 - 2A 2 , (78) 

and the equivalence of the cross-derivative condition (Y x ) t — (Y t ) x to the 
mKdV equation ([38|) for W proves that one has obtained a Darboux involution 
and a Lax pair, with the correspondence y = BY. 

The auto-BT of mKdV is obtained by the elimination of Y, i.e. by the 
substitution 

Log(BY) = a' 1 J(w- W)dx (79) 



in the two equations (|72[)-(73) for the gradient of y — BY. 

The singular manifold equation, defined (9^] as the constraint put on (p for 
the truncation to exist, is obtained by the elimination of W between S and 
C, 

bC-S-6\ 2 = 0, (80) 

and it is identical to that of the KdV equation. 

Remark. The fact that, in the Laurent series (|60|), mt(0) (the "l.h.s. ") 
and ut(oo) (the "constant level coefficient") are both solutions of the same 
PDE is not sufficient to define a BT, since any nonintegrable PDE also enjoys 
this feature. It is necessary to exhibit both the Darboux involution and a 
good Lax pair. 
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Most 1 + 1-dimensional PDEs with the PP have been successfully pro- 
cessed by the singular manifold method, including the not so easy Kaup- 
Kupershmidt |t3] and Tzitzeica |p8| equations. 

The extension to 2 + 1-dimcnsional PDEs with the PP has also been in- 
vestigated (Ref. (i^] and references therein). 

4.3 Singlevalued solutions of the Bianchi IX cosmologi- 
cal model 

Sometimes, the no-log conditions generated by the test provide some global 
information, which can then be used to integrate. 

The Bianchi IX cosmological model is a 6-dim system of three second order 
ODEs 

(Log A)" =A 2 -(B- Cf and cyclically, ' = d/dr, (81) 
or equivalently 

(Logwi)" = + — ^2^3/^1' A = &2Va/wi, u\ = BC and cyclically. 

(82) 

One of the families ||, |lj 

A = x^ 1 + o-2X + 0(x 3 ), X = t-t 2 , 

B = X - 1 +b 2 X + 0(x 3 ), (83) 

C = X - 1 +C2X + 0( X 3 ), 

has the Fuchs indices —1,-1,-1,2,2,2, and the Gambier test detects no 
logarithms at the triple index 2. The Fuchsian perturbative method 

N 2+JV— n 

A = X - 1 J2 £H E af X \x = r-r 2 , and cyclically, (84) 

n—0 j — — n 

then detects movable logarithms at (n, j) — (3, —1) and (5, —1) |pl[ , and the 
enforcement of these no-log conditions generates the three solutions : 

(4 0) = c 2 0) and b^l = c^l) or cyclically, (85) 

4°) = 4°) = 4°) = 0, (86) 

a^l = b^l = c% (87) 

These are constraints which reduce the number of arbitrary coefficients to, 
respectively, four, three and four, thus defining particular solutions which 
may have no movable critical points. 

The first constraint (85) implies the equality of two of the components 
(A, B, C), and thus defines the 4-dim subsystem B = C |h|, whose general 
solution is single valued, 

A _ fci B-C- fc i sinnfc i( T - T i) 



sinh ki (t - T\ ) ' k\ sinh fc 2 (r — T2 ) 
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The second constraint (|8g) amounts to suppress the triple Fuchs index 
2, thus defining a 3-dim subsystem with a triple Fuchs index —1. One can 
indeed check that the perturbed Laurent series J8^) is identical to that of the 
Darboux-Halphen system f37| , |52|| 

u>[ = u>2Ll>3 — UJ1W2 — ^1^3, and cyclically, (89) 

whose general solution is single valued. 



The third and last constraint (87) amounts to suppress two of the three 



Fuchs indices —1, thus defining a 4-dim subsystem whose explicit writing is 
yet unknown. With the additional constraint 

a<°>+&< 0) +c< 0) =0, (90) 

the Laurent series ( |S3|) is identical to that of the 3-dim Euler system (1750) 
Q| , describing the motion of a rigid body around its center of mass 

lu[ = and cyclically, (91) 

whose general solution is elliptic (J] 

ujj = (Log(p(r - r , 52, 33) - e,))', j = 1,2,3, (t , 32,53) arbitrary, (92) 

p' =4(p-ei)(p-e 2 )(p-e 3 ) = 4p 3 -. 9 2p-33- (93) 

The 4-dim subsystem (the one without (|90|)) defines an extrapolation to four 
parameters of this elliptic solution, quite probably single valued, whose closed 
form is still unknown. 

One thus retrieves by the analysis all the results of the geometric assump- 



tion of self-duality 1 48 , even slightly more 



4.4 Polynomial first integrals of a dynamical system 

A first integral of an ODE is by definition a function of x,u(x),u'(x), . . . 
which takes a constant value at any x, including the movable singularities of 
u. Consider for instance the Lorenz model 

dx . . dy dz . . . . 

— — = a\y — x), — — = rx — y — xz, —— = xy — bzix — y). (94) 
dt vy h dt y ' dt y v y ' y ' 

First integrals in the class P(x, y, z)e xt , with P polynomial and A constant, 
should not be searched for with the assumption P the most general polynomial 
in three variables. Indeed, P must have no movable singularities. The movable 
singularities of (x, y, z) are 

x~2ix _1 , V 2icr~ 1 x~ 2 , z 2ax~ 2 , indices (-1, 2, 4), (95) 

therefore the generating function of admissible polynomials P is built from 
the singularity degrees of (x, y, z) |62| 

1 



(1 — ax)(l — a 2 y){l — a 2 z) 



1 + ax + a 2 (x 2 + y + z) + a 3 (x 3 + xy + xz) 



a 



(x 4 + x 2 y + x 2 z + yz + z 2 + y 2 ) + . . . (96) 
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defining the basis, ordered by singularity degrees, 

(1), (x), (x 2 ,y,z), (x 3 ,xy,xz), (x 4 ,x 2 y,x 2 z,yz, z 2 ,y 2 ), ... (97) 

The candidate of lowest degree is a linear combination of (x 2 ,y,z), which 
indeed provides a first integral 



K x = {x 2 - 2az)e 2at , b = 2a. 



(98) 



Six polynomial first integrals are known [^0 with a singularity degree at most 
equal to four, and these are the only ones [63|. 

4.5 Solitary waves from truncations 

If the PDE is nonintegrable or if one only wants to find particular solutions, 
the singular manifold method of Section 4.2 still applies, it simply produces 
less results. For autonomous partially integrable PDEs, the typical output 
is a set of constant values for the unknowns S,C, A, B, R in the determining 
equations (|77j). In such a case, quite generic for nonintegrable equations, the 
integration of the Riccati system (^2|) yields the value 



X 



-tanh -(£-&), £ = x - ct, k 2 = -25, c = C, 



(99) 



the singular part operator T> has constant coefficients, therefore the solutions 
r in ( |7l| ) are solitary waves r = /(£), in which / is a polynomial in sechfc£ 
and tanh fc£. This follows immediately from the two elementary identities EG] 



tanh z — 



1 



tanh ; 



-2i sech 



2z + i- 



, tanh z 



1 



tanh ; 



2 tanh 



2z- 



(100) 

In the simpler case of a one- family PDE, the (degenerate) Darboux invo- 
lution is 



u = V Log i\i + U, d x Log -0 = X 1 j 



(101) 



and the above class of solitary waves r = /(£) reduces to the class of poly- 
nomials in tanh(fc/2)£. In the example of the chaotic Kuramoto-Sivashinsky 
(KS) equation 



u t + uu x + im xx + bu 

xxx ~r VUxxxx — 

0, i/^0, 



one finds 



15(16^- 6 2 ) 
76^ 



d x , 



u = 2? Log cosh -(£-£ ) + c, (c,^o) arbitrary, 



(102) 

(103) 
(104) 



in which b 2 l{p.v) only takes the value s 0, 144/47, 256/73, 16, and k is not ar- 
bitrary. In the quite simple writing (|l04j) , much more elegant than a third 
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degree polynomial in tanh, the only nonlinear item is the logarithm, T> be- 
ing linear and cosh solution of a linear system. This displays the enormous 
advantage to take into account the singularity structure when searching for 
such solitary waves. 

See e.g. Ref. |25) for a recent application to coupled Ginzburg-Landau 
equations. 

4.6 First degree birational transformations of Painleve 
equations 



At first glance, it seems that the truncation procedure described in section 4.2 
should be even easier when the PDE reduces to an ODE. This is not the case, 
because, in addition to the Riccati variable % or y of the truncation, there 
exists a second natural Riccati variable and therefore a homographic depen- 
dence between the two Riccati variables, which must be taken into account 
under penalty of failure of the truncation. 

Indeed, any Nth order, first degree ODE with the Painleve property is 
necessarily |7^, pp. 396-409] a Riccati equation for U^ N ~ 1 \ with coefficients 
depending on x and the lower derivatives of U, e.g. in the case of P6, 

U" = A 2 (U, x)U 12 + Ai{U, x)U' + A {U, x, A, B, T, A). (105) 

Then the Riccati variable of the truncation (denote it Z) is linked to U' by 
some homography, 

(U' + g 2 )(Z- 1 - g x ) - g = 0, g ? 0, (106) 

in which 30,31,32 are functions of (U, x) to be found. Implementing this 
dependence in the truncation |2^ ] provides a unique solution for P6, which 
is the unique first degree birational transformation, first found by Okamoto 
|7l, 



N = x(x-l)U' 6 9i e x -l 
-U U(U-l)(U-x) U U-l U-x 

x(x - i)u' 6 e 1 e x -i 



(107) 

1 ■ (108) 
u(u—\)(u — x) u u — 1 u — x 

Vj = oo,0,l, x : (6) +6| - (N/2) 2 ) 2 - (20j0j) 2 = 0, (109) 

n = ^(^_e 2 fe ), (no) 

with the classical definition for the monodromy exponents, 

e 2 x = 2a, e 2 = -2p, e 2 = 2 1 ,e 2 x = \- 25, (111) 

= 2A, e 2 = -2B, e 2 = 2r, e 2 = 1 - 2A. (112) 

The equivalent affine representation of ( |l09] )-( |iT(]| ) is 

0j = i-5(E ^) + 5' i = ^-5(E^) + 5' ( 113 ) 

n = i-^e fe = -i + ^0 fc = 2(^-e J -), j = 00, 0,1, x, (114) 
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in which j, k = oo, 0, 1, x. 

The well known confluence from P6 down to P2 then allows us to recover 
[2^1 all the first degree birational transformations of the five Painleve equa- 
tions (PI admits no such transformation because it depends on no parameter), 
thus providing a unified picture of these transformations. 

5 Liouville integrability and Painleve integra- 
bility 

A Hamiltonian system with N degrees of freedom is said Liouville-integrable 
if it possesses N functionally independent invariants in involution. In gen- 
eral, there is no correlation between Liouville-integrability and the Painleve 
property, as seen on the two examples with N = 1 

H(q,p,t) = ^--2q 3 -tq, (115) 
2 

H(q,p,t) = ^-+q 5 , (116) 

in which the first system is Painleve-integrable and not Liouville-integrable, 
and vice versa for the second system. However, given a Liouville-integrable 
Hamiltonian system, which in addition passes the Painleve test, one must try 
to prove its Painleve integrability by explicitly integrating. 
Such an example is the cubic Hcnon-Heiles system 

H = ^{pi+p 2 2+c 1 q^+c 2 ql)+aq 1 ql-^f3ql + -c 3 q^ 2 , a ^0,(117) 



q'{ + c iqi - pq{ + aql = 0, (118) 

q-2 + c 2 q 2 + 2a qi q 2 - c 3 q- 3 = 0, (119) 

which passes the Painleve test in three cases only, 

(SK) : 0/a = -l,CL = C2, (120) 

(K5) : 0/a = -6, (121) 

(KK) : 0/a = -16, ci = 16c 2 . (122) 



In these three cases, the general solution q\ (hence g|) IS indeed singlevalued 
and expressed with genus two hyperelliptic functions. This was proven by 
Drach in 1919 for the second case, associated to KdV 5 , and only recently Q 
in the two other cases. This proof completes the result of Ref. Q, who found 
the separating variables (a global object) by just considering the Laurent 
series (a local object), following a powerful method due to van Moerbeke and 
Vanhaecke ||. 
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6 Discretization and discrete Painleve equa- 
tions 



This quite important subject (the integrability of difference equations) is re- 
viewed elsewhere in this volume p9[ , so we will just write a few lines about 
it, for completeness. 

Let us consider the difference equations or ^-difference equations (we skip 
for shortness the elliptic stepsize 



Vx Vft : E(x, h, {u(x + kh), k - k = 0, . . . , N}) = 0, (123) 
VxVg : E(x,q,{u(xq k ), k - k = 0, . . . , N}) = 0, (124) 

algebraic in the values of the field variable, with coefficients analytic in x and 
the stepsize h or q. As compared to the continuous case, the main missing item 
is an undisputed definition for the discrete Painleve property. The currently 
proposed definitions are 

1. E3] There exists a neighborhood of h = (resp. q = 1) at every point 
of which the general solution x — *■ u(x,h) (resp. x — > u(x,q)) has no 
movable critical singularities. 

2. [j] The Nevanlinna order of growth of the solutions at infinity is finite. 

but none is satisfactory. Indeed, the first one says nothing about discrete 
equations without continuum limit, and the second one excludes the continu- 
ous P6 equation. 

Despite the lack of consensus on this definition, a discrete Painleve test 
has been developed to generate necessary conditions for the above properties. 
Of exceptional importance at this point is the singularity confinement method 
|51{ |, which tests with great efficiency a property not yet rigorously defined, 
but which for sure will be an important part of the good definition of the 
discrete Painleve property. The approach developed by Ruijsenaars [ |B5| for 
linear discrete equations, namely to require as much analyticity as possible, 
should be interesting to transpose to nonlinear discrete equations. 

Just for consistency, an interesting development would be to display a 
discrete version of (|T^) escaping all the methods of the discrete test. 

Let us say a word on the discrete analogue of the Painleve and Gambier 
classification. These second order first degree continuous equations all have 
a precise form (u" is a second degree polynomial in u' , the coefficient of u' 2 
is the sum of at most four simple poles in u, etc), directly inherited from 
the property of the elliptic equations isolated by Briot and Bouquet. In the 
discrete counterpart, the main feature is the existence of an addition formula 
for the elliptic function p of Weierstrass. As remarked earlier by Baxter and 
Potts (see references in f2jj| ), this formula defines an exact discretization of (Q). 
Then, all the autonomous discrete second order first degree equations with the 
(undefined!) discrete PP have a precise form resulting from the most general 
discrete differentiation of the addition formula, and the nonautonomous ones 
simply inherit variable coefficients as in the continuous case. Of course, the 
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second order higher degree (mostly multi-component) equations are much 
richer, see details in the review ppf . 

Another open question concerns the continuum limit of the contiguity re- 
lation of the ODEs which admit such a relation. The contiguity relation of 
the (linear) hypergeometric equation has a continuum limit which is not the 
hypergeometric equation, but a confluent one. On the contrary, the contiguity 
relation of the (linearizable) Ermakov equation has a continuum limit which is 
again an Ermakov equation |5j| . One could argue that the latter depends on 
a function, and the former only on a finite number of constants. Nevertheless, 
this could leave the hope to upgrade from P5 to P6 the highest continuum 
limit for the contiguity relation of P6 |7|, ITE, [2(| . 



7 Conclusion 

The allowed space forced us to skip quite interesting developments, such as 
the relation with differential geometry j^] , or the way to obtain the nonlinear 
superposition formula from singularities |73| ], or the weak Painleve property 
(7§ Lecons 5-10,13,19] g |§. 

For applications to nonintegrable equations, not covered in this Special 
issue, the reader can refer to tutorial presentations such that |l8|, |67|. 
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